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1 Introduction 


The quark model predicts the existence of heavy baryons composed of single, double and 
triple quarks. Essential improvement has been achieved in heavy baryon spectroscopy in 
resent years. All baryons with a single charm quark that are predicted by the quark model 
have been observed in experiments. Moreover, heavy baryons with a single bottom quark, 
such as Af,, Sfo, Ei, and Qb have also been discovered (for a review see [1]). 

These progresses in experiments stimulated future investigation for the properties of 
these baryons at LHC, as well as further theoretical studies on this subject. 

Remarkable information about the internal structure of baryons can be gained by study¬ 
ing their magnetic moments. Magnetic moments of the heavy baryons have long been under 
the focus of theoretical physicists, and they have been studied in framework of various ap¬ 
proaches. So far, magnetic moments of charmed heavy baryons have been calculated in 
framework of the naive quark model [2,3], relativistic quark model [4], chiral perturba¬ 
tion model [5], hypercentral model [6], soliton model [7], skyrmion model [8], bag model [9], 
QCD sum rules method [10], nonrelativistic quark model [11], phenomenological relativistic 
model [12], quark model with confinement law potential [13], and chiral constituent quark 
model, respectively. 

Magnetic moments of heavy hadrons have also been studied in numerous works within 
the QCD sum rules method. Magnetic moments of the A^ and baryons have been 
calculated in framework of the traditional QCD sum rules method in [10]; of the Ag and 
Sq {Q = c or b) baryons within the light version of the QCD sum rules method in [15,16], 
respectively. It should be noted here that magnetic moments of the spin-3/2 heavy baryons 
have already been analyzed within the same approach in [17]. 

The goal of the present work is calculation of the magnetic moments of heavy Eg, Sg 
and Dg baryons within the light cone QCD sum rules method. 

The paper is organized as follows. In section 2, sum rules for the magnetic moments 
of the above-mentioned baryons are constructed. Numerical results of our calculations are 
presented in section 3. This section further contains comparison of our results with the 
predictions of other approaches, and concluding remarks. 


2 Light cone sum rules for the sextet heavy baryons 

We start this section by giving a brief summarizing of the classihcation of heavy baryons 
in SU(3) symmetry. According to SU(3) symmetry, heavy baryons with single heavy quark 
belong to either symmetric sextet or antisymmetric anti-triplet flavor representations. As 
has already been noted , magnetic moments of the Ag and Sg baryons belonging to anti¬ 
symmetric flavor representations of SU(3) are calculated within the light QCD sum rules 
method in [15] and [16], respectively. In the present work we calculate magnetic moments 
of the heavy baryons belonging to sextet representation of SU(3) group. 

In order to calculate magnetic moments of the heavy sextet baryons, we start by con¬ 
sidering the following correlation function: 

n(p,g)=i j |T{r/g(a;)r/g(0)}|0) , (1) 
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where rjg is the interpolating current of the heavy spin-1/2 baryon. The main task in 
constructing the sum rules for magnetic moments of the heavy sextet baryons is the calcu¬ 
lation of the correlation function in terms of the photon distribution amplitude by using the 
operator product expansion (OPE) over twist from one side, and in terms of the hadrons 
from the other side, and then equating both representations. Calculation of the correlation 
function from the hadronic side is accomplished by inserting a complete set of hadrons car¬ 
rying the same quantum numbers of the interpolating current rjq. Isolating the the ground 
state’s contribution, we obtain. 


n 


(0|hQ|-BQ(p2)) 

pI - 


{B{p2h{q)\BQ{pi))^ 


{Bq{pi)\vq\0) 

pI - 


+ • • • , 


( 2 ) 


where dots correspond to the contributions of higher states and continuum, p\ = P 2 + <?• In 
further analysis we shall make the replacement p 2 = p. The matrix elements in Eq. (2) are 
determined as. 


{'^\Vq\Bq{p)) = \qUbq{p) , 


{Bq{p)\Bq{pi))^ = 


/i7m - /2 


za 


fJ-i' 


2m 


Bq 


q 


UBniP 


•gm) > 


(3) 

(4) 


where is the photon polarization vector, fi and /2 are the form factors. Using the 
equation of motion, Eq. (4) can be written as. 


{Bq{p)\Bq{pi))^ 


ubq{p) 


(/l + /2)7 m + 


(p + Pi)m 

2mBQ 


UBQ{pi)e^ ■ 


(5) 


Since the photon involved in these transitions is a real photon , we set = 0, and hence 
in calculation of the magnetic moments of the heavy sextet baryons the values of the form 
factors are needed only at = 0. 

Substituting Eqs. (3) and (5) in Eq. (2), and performing summations over spins of the 
heavy baryons, for the hadronic of the correlation function we get. 


n 


3^2 + 

Bq p2 _ 


(/l + + 


{P> + Pl)l^ 
2mBQ 


ih + ^Bg) 

pI - bt-Iq 


( 6 ) 


We observe from Eq. (6) that the correlation function contains many structures, any of 
them can be chosen in calculating magnetic moments of the sextet baryons, and in this 
respect we choose the structure This structure envelopes the magnetic form factor 
fi + / 2 , and at g^ = 0 it gives the magnetic moment pbq of the heavy baryons in units 
of eh/2mBQ- As a result, the correlation function can be written in terms of the magnetic 
moment of heavy baryons as, 

2 1 1 

-{p + qT 


(7) 


In order to calculate the correlation function in terms of quark and gluon degrees of 
freedom and photon distribution amplitudes, the expressions of the interpolating currents 
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of the heavy baryons are needed. The general form of the interpolating currents of the 
heavy spin-1/2 positive parity baryons is given as (see for example [18]), 

where a, b, c are the color indices, C is the charge conjugation operator, and t is an arbitrary 
parameter whose value at t = — 1 gives the so-called Ioffe current. The quark contents of 
the sextet heavy sextet baryons are given in Table 1. 
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Table 1: Quark contents of the heavy sextet baryons. 

Using the expression for the interpolating current and Wick’s theorem, the theoretical 
part of the correlation function for the can be written as, 

= -3|755'd(a;)5'b(a;)S'„(a;)75 -f 'y5Su{x)Sl{x)Sd{x)'y5 + 'y^Su{x)'^^Tr 

+ 75*S'd(a;)75Tr [S'„(a;)S'b(a;)] + t(^^Sd{x)'^^S'^,{x)Su{x) + '^^Su{x)'-i^S'^{x)Sd{x) 

+ Sd{x)S'^{x)'^^Su{x)'-f^ + Su{x)S'^,{x)'-i^Sd{x)'^^ + '-f^Su{x)Tr [Sb{x)'y5Sd{x)] 

+ Su{x)'-i^Tr [Sb{x)Sd{x)'y5] + 'y^Sd{x)Tr [S'«(a;)755'b(a;)] Sd{x)'-ibTr [5'«(a;)S'b(a;)75] j 

-h t^{Sd{x)'^bS'b{x)'ybSu{x) + Su{x)'-i^S'^{x)'-ibSd{x) + Sd{x)Tr [5'fe(a;)75S'/(a;)75] 

-h S'„(a;)Tr [S'd(a;)755'fe(a;)75] , 

where S' = CS'^C, T symbolizes transposition operator, and S is the quark (light or 
heavy) propagator. The corresponding expressions of the correlation functions for the 
other members of the sextet baryons can be found from the by making the following 
replacements: 

ns? = (m ^ d), 

ns++ = (d ^ m) , 
nw = d ), 

= u^^d^u) , 

jf.'+ ^ (d ^ s), 

11“=° = {u ^ d,d ^ s) , 
n“6° = s), 

n = n ** (ti — y d, d — y s') , 

{u^s,d^s), 
nw _ . 
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The correlation function described by Eq. (9) contains three different parts: a) Pertur¬ 
bative part, when a photon is radiated from short distance (perturbative contribution); b) 
a photon is radiated from short distance from the quark propagators, and light quarks form 
a condensate; c) a photon is radiated from long distance (nonperturbative contribution). 

In order to calculate the perturbative contribution of the correlation function, it is 
enough to make the replacement 

J - y) , 

for one of the quark propagator, where Sfree is the free quark operator; and the other 
two are being the free propagators. The expressions of the free light and heavy quarks in 
coordinate representation are given as, 

fUq 

47r2 * 47r2 (7^)2 

where Ki are the Bessel functions. 

The contribution of part (c) can easily be calculated by replacing one of the light quark 
operator with. 


^free _ 


gfree _ 

Q 


where Tj are the full set of Dirac matrices, and other quark operators that involve pertur¬ 
bative, as well as nonperturbative terms. The explicit forms of the “full” quark propagators 
can be found in [19] and [20]. 

Nonperturbative contribution is realized as the matrix element of the nonlocal operator 
qViQ between the vacuum and one-photon states.These matrix elements are described in 
terms of photon distribution amplitudes, as are given below (see [21]): 


{7(g)|g(a;)a^^g(0)|0) = -ieqqq{e^q^ - e^q^) / ( x^-yiu) + 


X 


2{qx 


-eq{qq) 


ex 


Xu [Sii- q^. 


qx 


ex\ 

x^,\eu- qu — 
qx J 


due^^^^h^iu) 


{j{q)\q{x)j^q{0)\0) = Cqfs^ (sf, - q^,—^ [ 

V QX/ Jq 

(7(g)|g(a;)7/,75g(0)|0) = -^Cqf^^equays''q""x^ 

{l{q)\q{x)gsGf,u{vx)q{0)\0) =-ieq{qq) {Ef^qu - Euqi,) J 

{xiq)\qix)gsGf,J'y5{vx)q{0)\0) = -ieq{qq) {ef,q„ - e^q,,) [ 
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Sfji - q^l—] (gau - —{qaXu + qyXa) 1 ?/? 

qx J \ qx ' 


- e 


(7(g)|g(a;)5(^G/,^(ra)7c,75g(0)|0) = eqf2,^q^{e^q^ - e^q^) j 
{'l{q)\q{x)gsG^u{vx)i'^aq{^)\^) = egf^^q^ie^q^ - e^qf,) [ 
{l{q)\q{x)(Tai3gsGf,^{vx)q{0)\0) = eq{qq) 

- 9m— ) (gpu- — [q^Xy + qyXp) ) g« 

qxj\ qx J 

- 9a^^ {gup - + 9/32/,.)^ 9m 

~ ~ 9aa;;,)^ (?,. 


Vaie^G,+^^<,)q^q-i{ai) 


+ [ep- qp—] (gua - —{quXa + qaXu) ) q^, 
' qx J \ qx ' 


Vaie^G,+v^,)<i-j-^^ai) 


H- {qi^Xu - quX^,){eaqp - £pqa) / 

qx J 

+ ~ qpXa){£fMqu - £uq,,) j PQ;ie*("«+^“»)''"^71(ai) | . 

In the definitions given above, x is the magnetic susceptibility of the quarks, is the 

leading twist-2, 'ip'"{u), 'ip°‘{u), A and V are the twist-3, and h^{u), A, % {i = 1, 2, 3, 4) 
are the twist-4 photon DAs, respectively, whose explicit expressions are given in Appendix 
A. The measure Vai is defined as 

/ Dai = / daq / daq / dag6{l — aq — aq — ag). 

Jo Jo Jo 

In constructing the sum rules for the magnetic moments of heavy sextet baryons it is 
necessary to equate the coefficients of the structure to from the phenomenological and 
theoretical representations of the correlation function. The following steps in obtaining the 
final form of the sum rules for the magnetic moment are: Fourier transformation, Borel 
transformation over the variables and {p + qY variables, and continuum subtraction 
in order to suppress the contribution of the higher states and continuum. After these 
procedures we get the sum rules for the magnetic moment of the sextet heavy baryons, 
which can schematically be written in the following form, 

\2 ,, _ ns,tw 

-n ’ 

where the function j{E,theor ^Q^tain perturbative and nonperturbative contributions. 
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The expression for Yi^^iheor jg quite lengthy, and for this reason we do not present it 
here. 

In order to determine the magnetic moments of heavy sextet baryons, the value of 
the overlap amplitude is needed, which can be found from the analysis of the two-point 
correlation function. The residues of the heavy sextet baryons are calculated in [22], and 
their expressions can be found in this work. 

3 Results and discussion 

Having obtained the sum rules for the magnetic moments of the heavy sextet spin-1/2 
baryons, we are now ready to perform the numerical analysis. The main input param¬ 
eters of the light cone QCD sum rules for the magnetic moments are the photon distri¬ 
bution functions DAs [21]. Their explicit expressions are presented in Appendix A. The 
values of the other input parameters are: {uu) ^=iGeV = {dd)^=iGeV = —(0.243)^ GeV^, 
{ss)^,=iGeV = 0.8{uu)^=,GeV [23], mg = (0.8 ± 0.2) [22], A = (0.5 ± 0.1) GeH [24] 

/ 3 ^ = -0.0039 GeH2 [21], m,(^= 2 GeV) = (HI ± 6) MeV [25]. 

Few words about the magnetic susceptibility involved in the numerical analysis are in 
order. Its value is estimated to have the value x(l GeV) = —4.4 GeV~^ in [26]. In [21], 
using the vector dominance model ansatz and QCD sum rules its value is predicted to be 
x(l GeV) = —(3.15±0.15) GeV~^. Furthermore, from an analysis of the radiative decays of 
heavy mesons its value is found to be x(l GeV) = —(2.85 ±0.50) GeV~^ [27]. In numerical 
analysis we have used all these predicted values of the magnetic susceptibility. 

It should also be remembered that, the sum rules for the magnetic moments involve 
the following auxiliary parameters: The arbitrary parameter t appearing in the expressions 
of the interpolating currents, Borel mass parameter M^, and the continuum threshold Sq. 
The magnetic moment must be independent of these parameters. In order to hud “regions” 
of these parameters where magnetic moment exhibits good stability with respect to their 
variations we proceed as follows. Firstly, we try to find the upper and lower bounds of Borel 
mass parameter at hxed values of Sq and t. The upper bound of can be determined 
by requiring that the contribution due to the continuum threshold should be less than 
half of the contribution coming from the perturbative part. The lower bound of can 
be found by demanding that the higher powers of 1/M^ contribute less than the leading 
twist contributions. Our numerical analysis shows that these two conditions are satished 
simultaneously if ranges in the following regions: 

2.0 GeI/2 < < 3.0 GeH^ 

2.2 GeV^ < < 3.4 GeH^ 

2.5 GeH^ < < 4.0 GeH^ 

5.0 GeV^ <M^ < 6.0 GeV^ Sb 

5.0 GeV^ <M^ < 6.4 GeV^ 

5.2 GeV^ <M^ < 7.0 GeV^ Db 

In regard to the continuum threshold Sq appearing in the sum rules, it is known that 
the values of this arbitrary parameter is related to the energy of hrst excited state. The 
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difference ^/^—mground, where rUground is the ground state mass of the baryon, is the energy 
needed to excite the particle to its hrst energy state. This quantity usually varies between 
0.3 GeV and 0.8 GeV. Analysis of the mass sum rules shows that, in order to reproduce 
the experimental mass of the sextet baryons the continuum threshold should vary in the 
following domain: 

(3.1 ±0.1) Gel/ Sc 
(3.2 ±0.1) Gel/ S'c 
(3.3 ±0.1) Gel/ fic 
(6.6 ± 0.2) Gel/ Sf, 

(6.7 ± 0.2) GeV 
(6.8 ± 0.2) GeV 

Having determined the sum rules, input parameters, and working regions of all auxiliary 
parameters, we perform numerical analysis to calculate the magnetic moment ^Bq of the 
heavy sextet spin-1/2 baryons, whose results we can summarize as follows. As examples, 
in Figs. 1 and 2 we present the dependence of /isg on at several hxed values of t 
and at the hxed value of sq listed above for the baryons S^"'' and S°, respectively. We 
observe from these hgures that the magnetic moments of Sq show good stability with 
respect to the variation of the Borel mass parameter in aforementioned domains. The next 
step in hnding the values of the magnetic moments of the baryons under consideration is to 
determine working region of the arbitrary parameter t, where ^Bq is practically independent 
of its variation. For this purpose, in Figs. 3 and 4 the dependence of /isg on cos 6*, 
where tan 6* = t for the and baryons, at several hxed values of the Borel mass 
parameter and at the predetermined value of the continuum threshold sq, are presented, 
respectively. It follows from these hgures that the magnetic moments of the heavy sextet 
spin-1/2 baryons seem to be practically independent of the arbitrary parameter t in the 
domain —0.9 < cos 9 < —0.6, and insensitive to the the chosen values of Sq. Our hnal results 
on the magnetic moments of the heavy sextet spin-1/2 baryons are presented in Table 2. 
For comparison, in the same Table we present the predictions of other approaches on the 
magnetic moments of the relevant baryons, such as nonrelativistic quark model [11], bag 
model [9], phenomenological relativistic quark model [12], quark model with conhnement 
law potential [13], relativistic quark model [4], hypercentral model [6], chiral constituent 
model [14], and QCD sum rules method [10]. 

We observe from this table that almost all approaches give, more or less, similar pre¬ 
dictions, except the results of [9] and [4] (especially for the charmed baryons), which are 
smaller. 

Using our results one can easily deduce the following relations among the magnetic 
moments of heavy baryons. 

Fs++ + h-sS? — 2p-s+ ) 
h'S^+ ± h-og — , 

A 2/i=/y — /isg ± , 

Fe+ a Fs" — 2Ats° ) 

Fs+ A — 2p.=^o , 
hE+ A 2/i2'- — /is- A 2/i2'o , 
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[4] 

[6] 

[9] 

[10] 

[11] 

[12] 

[13] 

[14] 

Our Work 


-1.04 

-1.015 

-1.043 

-1.6 ± 

0.2 

-1.37 

-1.38 

-1.39 

-1.6 

-1.5 


0.36 

0.5 

0.318 

-H 

o 

0.1 

0.49 

0.49 

0.525 

0.3 

0.5 


1.76 

2.279 

1.679 

2.1 ± 

CO 

o 

2.35 

2.36 

2.44 

2.22 

2.4 


-0.95 

-0.966 

-0.914 

- 


-1.18 

-1.12 

-1.12 

-1.32 

-1.2 


0.47 

0.711 

0.591 

- 


0.89 

0.75 

0.796 

0.76 

0.8 

HO 

-0.85 

-0.96 

-0.774 

- 


-0.94 

-0.86 

-0.85 

-0.9 

-0.9 


-1.01 

-1.047 

-0.778 

- 


-1.22 

- 

-1.256 

- 

-1.3 


0.53 

0.592 

0.422 

- 


0.64 

- 

0.659 

- 

0.6 


2.07 

2.229 

1.622 



2.5 

- 

2.575 

- 

2.4 


-0.91 

-0.902 

-0.66 

- 


-1.02 

- 

-0.985 

- 

-1.2 

-/O 

0.66 

0.766 

0.556 

- 


0.9 

- 

0.930 

- 

0.7 


-0.82 

-0.96 

-0.545 

- 


-0.79 

- 

-0.714 

- 

-0.8 


Table 2: Magnetic moments of the heavy sextet, = i"*" baryons in units of the nuclear 
magneton 


The direct measurement of the magnetic moments of heavy baryons are unlikely in the 
near future. Therefore, any indirect estimations of the magnetic moments of the heavy 
baryons could be very useful. For example, it could help extracting information about the 
mass spectrum of the heavy baryons. As is well known, experimentally measured mass 
difference is attributed to the hyperhne splitting. Moreover, the magnetic moments of 
quarks are proportional to the chromomagnetic moments, which determine the hyperhne 
splitting in the baryon spectrum. Following this reasoning, the magnetic moments of the 
Ac and A;, baryons are estimated in [28]. Hence, this approach could be an essential tool in 
estimating the magnetic moments of the heavy baryons. 

In conclusion, the magnetic moments of the heavy sextet baryons are calculated 

in framework of the light cone QCD sum rules method. Empirically, linearly independent 
relations among the magnetic moments of the sextet baryons are obtained. Comparison of 
our results with the predictions of other approaches is presented. 




Appendix A: Photon distribution amplitudes 

Explicit forms of the photon DAs [21], 

r - 

= Quu 1 + (p 2 {l^)C 2 {u — u) , 

= 3[3(2m - if - 1] + ^(15^7 - 5w;f)[3 - 30(2m - 1)^ + 35(2m - 1)^] , 

<(«) = [1 - (2m - 1)1 [5(2m - 1)2 - 1] ^ , 

A{ai) = SGOttgaga^ 1 + wf^{7ag - 3) , 

V{ai) = 540wf{ag - afagagO^g , 
hr^{u) = — 10(1 + 2k^)C2 {u — u) , 

A{u) = 4:0 u‘^u‘^{3k, — + 1 ) + 8 (^ 2 ^ — 3 C 2 )[m-u (2 + 13uu) + 2m^(10 — 15m + 6uf ln(M) 

+ 2#(10 - 15m + 6m 2) ln(M)] , 

120(3((^2 a C 2 )(*^g O'g'jO'qO'qCYg , 

7i(ai) = 3Dal{aq - aq)[{K - m+) + (Ci - Cl)(l - 2ag) + ( 2(3 - Aag)] , 

'TliiOif 120(3((^2 C 2 )(*^9 Oiq)(y.q(y.q(y.g , 

Tfai) = 30q;2(q;- — aq)[{K + kA) + (Ci + Ci')(l — 20;^) + (^ 2(3 — ^oig)] , 

S{ai) = 30Q;2{(fi; + m''')(1 — otg) + (Ci + Ci~)(l — 0'g){l — 2ag) 

A C2[3(cig ^q) Clg(l *^ 9 )]} ) 

§{ai) = -30a2{(K - K+)(l - ag) + (Ci - Cl)(l - agfl - 2ag) 

A C2[3(ciq ^q) *^ 3(1 *^ 3 )]}' 

The parameters entering the above DA’s are borrowed from [21] whose values are 992(1 GeV) = 
0, wf = 3.8 ± 1.8, = -2.1 ± 1.0, M = 0.2, m+ = 0, Ci = 0.4, C 2 = 0.3, (1 = 0, and 

C 2 + = 0. 
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Figure captions 

Fig. (1) The dependence of the magnetic moment of the baryon on at several 
hxed values of t, and at sq = 9.0 GeV^. 

Fig. (2) The same as Fig. (1), but for the baryon, at sq = 41.0 GeV"^. 

Fig. (3) The dependence of the magnetic moment of the baryon on cos 6, at several 
hxed values of and at sq = 9.0 GeV^. 

Fig. (4) The same as Fig. (3), but for the baryon, at Sq = 41.0 GeV"^. 
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